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Abstract 

We study bound states embedded into the continuum of edge states in two-dimensional topological insulators. These 
states emerge in the presence of a short-range potential of a structural defect coupled to the boundary. In this case the 
edge states flow around the defect and have two resonances in the local density of states. The bound state in continuum 
(BIC) arises due to an interference of the resonances when they are close to the degeneracy. We find the condition 
under which the BIC appears, study the spacial distribution of the electron density, and show that the BIC has a helical 
structure with an electron current circulating around the defect. 
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1. Introduction 

The existence of bound states embedded into the con¬ 
tinuum of quantum states was proposed at the dawn of 
quantum mechanics by von Neumann and Wigner [T] for 
certain attractive potentials of very specific form. The 
formation of the bound state in continuum (BIC) can be 
considered as a result of a destructive interference of par¬ 
tial waves in the potential due to which the wave function 
vanishes at large distance. It is essential that the BICs 
are true eigenstates of the Hamiltonian and therefore are 
orthogonal to the eigenstates of the continuous spectrum. 
Because of a rather complicated form of the potential these 
states were regarded a long time as mathematical curiosi¬ 
ties. However, further developments have produced a better 
understanding of the kind of potential that can create such 
bound states. Friedrich and Wintgen [5] have pointed out 
that the resonant states play an important role in the BIC 
formation and proved that a BIC could emerge as two cou¬ 
pled resonant states connected to a continuum are driven 
into degeneracy by changing a parameter governing the 
resonant states. 

The advances in the nanostructure technology have pro¬ 
vided new possibilities to create the required potential and 
the interest to BICs has strongly increased due to their 
manifestations in quantum transport, optical phenomena, 
and potential applications for quantum information devices. 
Herrick [3] and Stillinger |4j have proposed a construction 
of the potential supporting BICs for epitaxial heterostruc¬ 
ture superlattices. The localized states with properties 
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very close to BICs have been realized by Capasso et al. [5] 
in semiconductor heterostructures. In recent years the 
BICs are widely studied also in photonic crystals and were 
observed experimentally [5]. 

Theoretical studies of the occurrence of BICs in low¬ 
dimensional structures have been investigated for a variety 
of systems such as quantum Hall effect and narrow-wire 
circuits [7], systems of coupled quantum dots [HHIQ] , open 
quantum billiards m, trilayer graphene flakes connected 
to nanoribbon leads m- A distinct mechanism of the 
BIC formation was recently proposed for a chiral quantum 
system coupled to leads with a continuum energy band [13] . 

In the present paper we address the question of how the 
unique properties of topological insulators (TIs) affect the 
possibility of the BIC formation and come to the conclusion 
that at least in the case of two-dimensional (2D) TIs the 
emergence of BICs is facilitated by the presence of the 
topological order so that a BIC appears in a simple potential 
produced by a structural defect or impurity. 

The distinctive feature of the 2D TIs is the presence of 
helical edge states with massless Dirac spectrum in the 
gap of bulk states M- In the edge states the electrons 
move along the boundary and their spin is locked to the 
momentum because of strong spin-orbit interaction. In 
recent works [T51ITC] we have studied the electronic states 
induced by nonmagnetic defects and found that a defect 
with a short-range potential creates two bound states in 
contrast to topologically trivial insulators where only one 
bound state exists at such a defect. These states differ in 
both the pseudospin structure and the spatial distribution 
of the particle density. It is also interesting that they arise 
for any sign of the defect potential. The presence of two 
states is caused by the fact that in 2D TIs there are two 
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mechanisms of the bound-state formation. One mechanism 
is conventional. The bound state is formed by a potential 
attracting the quasiparticles of one of the bands. Another 
mechanism is specific for TIs. It is caused by the formation 
of an edge state circulating around the defect similarly to 
the edge states near the boundary. One can say that the 
defect effectively creates a boundary condition for the wave 
function. 

When the defect is located at a finite distance from 
the boundary, the bound state is coupled with the edge 
states giving rise to the formation of the edge states flowing 
around the defect. These states have a continuous spectrum 
with resonances in the local density of states near the defect 
at the energy close to the bound state energy m- 

In the present paper we study the true bound states, 
which form due to the interference of these resonances. 
These states have a square integrable wave function and 
the energy embedded into the continuum of the edge states 
flowing around the defect. Because of strong spin-orbit 
interaction these BICs have a helical structure. 


2. Model and equations determining the BIC 


Consider a 2D TI, which contains a defect located at 
a distance yd from the boundary. The defect creates a 
potential V(x,y — yd) localized in a small region. We use a 
layout in which the a:-axis coincides with the boundary and 
the y-axis is directed into the TI. The electronic states in 
the 2D TI are described by the model Hamiltonian proposed 
by Bernevig, Hughes, and Zhang (BHZ) for HgTe/CdTe 
quantum wells m- This model has the Sz symmetry 
since it does not take into account a possible spin-orbit 
interaction arising from the lack of the inversion symmetry. 
The corresponding generalization seems not to be essential 
at this stage and will be done elsewhere. 

Thus, the Hamiltonian of the system is 


H = 


(h{k) + V{r) 
\ 0 


0 ^ 

h*(-k) + V{r)J ’ 


where k is the momentum operator and 


h(k) 


/'M — {B+D)k^ A{kx + iky) \ 
A{kx-iky) -M+{B-D)k^) 


( 1 ) 


( 2 ) 


with M, A, B and D being the parameters of the BHZ 
model. The topological phase is realized when MB > 
0 [miiH]. In practically important cases of the HgTe/CdTe 
and InAs/GaSb/AlSb quantum wells, the parameters 
M, B,D < 0, and A > 0. The basis set of the wave 
functions is {|i?i f), t)i l-^^i i) i)} where \Ei f) 
and \Ei /,) are superpositions of the electron states of s- 
type and the light-hole states of p-type with spin up and 
spin down; |i7i f) and \Hi J,) are the heavy-hole p-type 
states with opposite spins. In what follows we will restrict 
ourselves by considering the symmetric model where D = 0. 

As boundary conditions for the wave function, we impose 
'I'(a; = 0,y) = 0 at the boundary and suppose that the 


wave function d'(a;, y) does not diverge at the infinity [x —>■ 
±oo, y —>■ oo). 

A simplification arising from the Sz symmetry of the 
BHZ model is that the Hamiltonian is block diagonal and 
therefore it is enough to consider only one spin component. 

The eigenfunctions of the Hamiltonian 0 are found 
using the Fourier transform with respect to x variable and 
the Laplace transform for y variable. This method allows 
one to find exactly the eigenfunctions in the case of the 
short-range potential where the interaction radius is small 
compared to the characteristic length scale of the wave 
function. The procedure of the solution was developed 
previously m, so we present here only final equations 
without technical details. 

For convenience we use the dimensionless variables: 

e = E/\M\, {x',y'} = {x,y}y^MjB^j^A/^/MB, , , 

v{x',y') = V{x,y)/\B\, b = yd^/M/B , 

and in what follows we shall omit the prime in the variables 
x,y. 

The Fourier-Laplace transform 'k(A:,p) of the eigen¬ 
functions is expressed through two variables: $(fc;u,6) 
and 'k(u,5), where d>(fc;u,5) is the Fourier transform of 
the y-derivative of 'I'(a::,y) at the boundary y = 0 and 
'^{v,b) is the value of the wave function at the defect 
d>(w, 6) = 'I'(0, &). The equation that determines dt(fc,p) 
reads as: 


[e-h{k,p)]'S{k,p) = az<^{k;v,b)+l 2 x 2 v{k,p)e ^P'i>{v,b). 

( 4 ) 

where 12x2 is the identity 2x2 matrix, az is the Pauli 
matrix, v{k,p) is the Fourier-Laplace transform of v{x,y). 
For simplicity, we will present hereafter the equations in a 
reduced form for v(k,p) « v. 

'^{v,b) is determined by the following equation: 

[I 2 x 2 - A(e; V, b)] «'(u, b) = C{e; u, b)E{e] b)x, (5) 


where lC{e]v,b) and E{e-,v,b) are 2x2 matrices which 
are expressed via the matrix elements of the Hamiltonian 
h(k); X is the following spinor: y = (1,-1)^; C(e;u, 6) is 
a function, which can be determined by the normalization 
of the wave function. 

The equations defining JC{e;v,b) and E{£]b) are rather 
cumbersome. In an explicit form they are given in|Appendix 

El 

The spinor ^(k;v,b) is determined by the following equa¬ 
tion: 


$(fc;u,6) = - ^ + C{£;v,b)xS[k - koie)], 

( 6 ) 

where ko{£) is the root of the determinant Ai(£, k) defined 
in Eq. (A.9|. It is easy to show from Eqs. (A.9), ( |A.6 ) and 
(A. III that fco(e) is exactly the spectrum of the edge states 
in the absence of the defect: 


fco(e) = -- 


( 7 ) 
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Eqs. ([^ and Q allow one to calculate the wave function 

oo c+ioo 

J 2 ^® J 27nA{s,k,p) 

— OO c—200 

X [Do{s,k,p)^{k-v,b) + v{k,p)e~^PDi{e,k,p)'i>{v,b)] . 

( 8 ) 

Let us turn to Eq. ([^ which plays a key role in finding the 
BIC. This equation has solutions of two kinds depending 
on the determinant of the matrix in the left-hand side 

A^(e;w, 6) = [l - lCii{e;v,b)] [l - IC 22 {e; v, b)] 

- ICi 2 (e;v,b)IC 2 i{e-,v,b). (9) 

If A^(e; V, b) ^ 0, the solutions of Eq. Q correspond to 
the continuum of the edge states flowing around the de¬ 
fect. If the determinant equals zero, this equation has one 
more solution which exists in the case where the function 
C{e-,v,b) is zero. This solution exists only at a discrete 
value of the energy eo defined by the root of the determi¬ 
nant: 

A^(eo;?^,&) = 0. (10) 

In this case 


'S{v,b) = Cbs 


1 


(l-ZCii) //Ci2 


( 11 ) 


e=eo 


with the constant Cbs being determined by the normaliza¬ 
tion. 

Eor such a solution to exist the wave function should be 
square integrable. Since ’^{x,y) is finite and tends to zero 
at y —oo, it is sufficient to require that 'k(x, j/) goes to 
zero when x —>■ ±oo faster than 

The asymptotics of 'I'(x,y) can be found from Eqs. (§ 
and It has the following form: 


dt)! —>■ oo, y) ~ 




8a 

OUe Qf, L 


P2 


-Do{e,k,-p2) 


o-piy 


Pi 


-Do{e,k,-pi) 


A'{e, k)B{e, k; b) 


k=ko(e) 


( 12 ) 


At this point we will take into account the following prop¬ 
erty of the matrices A'(e, k) and B{e, k\ b) in the integrand: 


A'{e,k)B{e,k;b){ J ) = 0, 


(13) 


which can be verified by direct calculation. Because of this 
equation we find that tl/(a;—>-oo,?/)—J-Oif 




(14) 


Combining Eq. (14) with Eqs. (10) and (11) we arrive 


at the following equation system that determines a bound 


state embedded into the continuum of the edge states 
flowing around the defect: 


1 - ICii{e;v,b) - ICi 2 {e-,v,b) =0 
1 - A22(e;u,6)-/C2i(e;i’,&) =0. 


(15) 


A solution of these equations, if it exists, defines the energy 
Sbs of the BIC and a relation between v and b under which 
this state exists. 


3. Energy of the BIC and conditions of its exis¬ 
tence 


In this section we show that Eq. (15) really has a solution 
and clarify how the energy of the BIC Sbs and the potential 
Vbs, at which it emerges, depend on the distance b for 
realistic systems. 


First we prove that Eq. (15) has one solution when the 


distance b is large enough. It is convenient to present the 


matrix /C(e; v, b) given by Eq. (A.l) in the form where the 
dependence on v is written explicitly 


/C(£; v,b) = V (s) -I- A(e; 6)], 


(16) 


Here the first term in the right-hand side is easily shown to 
be diagonal and the second term exponentially decreases 
with b. 

To begin we consider the limit where b ^ oo. In this 
case /Ci 2 = /C 21 = 0 and Eq. takes the form 


l-uJI?(^(e) =0 

1-vAf^ie) = 0 . 


(17) 


Here the first equation defines the energy of the electronlike 
bound states, Eeiy), and the second equation defines the 
energy of the holelike bound states, eh{v). The equations 
are compatible if ee{v) = eh{v). The energies Seiv) and 
£hiv) as functions of v were investigated in our previous 
works It was found that the equation £e{v) = eh(u) 

has a unique solution for both positive and negative v. 
Thus, the equation system © is satisfied when e = ±e 
and correspondingly v = ±h. 

When the distance b is finite but large enough so that 
^ 12,^21 ^ 1, one can seek a solution of Eq. (15) by 


expanding e and v near e and w at a given b. On putting 
£ = £ -|- i 5£ and v = v + Sv, we arrive at the following 
equations: 


-|- (5u = b) — Ai2(£; b)] 

^I?2<^£ -I- (5u = 'y^[Il22(c b) - J?2i(e; b )], 


(18) 


where A!i 2 = '5e-^ii,22(^) l^- H 1® ®^sy to show that the 
determinant of this system does not equal zero and there 
is a unique solution. 

A different situation occurs when the distance b is small. 


In this case the system (15) turns out to have no solution. 

Thus the BIC can exist if the distance b exceeds a thresh¬ 
old value and the potential v is equal to a definite value 
Vbs{b) for a given b. 
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Figure 1: (Color online.) Resonant energies as functions of the 
defect potential for a variety of distances between the defect and 
the boundary, b = 4,5,6, near the point where the resonances are 
degenerate. The asterisks show the energy of the BIC for a given b. 
The calculations were carried out for a = y/2. The inset illustrates 
the defect, the BIC and edge states flowing around the defect. 



Figure 2: (Color online.) The BIC energy (a) and the defect 
potential at which the BIC arises, (b) as functions of the distance 
between the defect and the boundary. The calculations were carried 
out for three values of the parameter a = 1.3, y/2^ 1.5. 


Now we present results of numerical studies of the BIC. 
Since the BIC arises due to the interference of the reso¬ 
nances in the local density of the edge states flowing around 
the defect, we present first the resonant energies in the 
vicinity of the avoided crossing point of the resonances 
originating from the electronlike and holelike bound states. 
Fig .[^ shows the resonant energies as functions of the defect 
potential v. The points where the BIC arises are shown by 
the asterisks for different b. 

The variation of the BIC energy with the distance be¬ 
tween the defect and the boundary is shown in Fig. [^a). 
Fig. [^b) presents the defect potential Vbs, at which the 
BIC arises, as a function of the distance b. It is seen that 
with increasing b the BIC energy and the defect potential 
tend to the limiting values £ and v defined by Eqs. (18). 

To complete the description of the BIC in 2D TI we 
present also the spacial distribution of the particle density 
in the BIC and the electron current density distribution. 
The wave function is given explicitly by Eq. ([^ where 
the matrices in the integrand should be calculated for 



Figure 3: (Color online.) (a) The spacial distribution of the particle 
density in the BIC induced by the defect located at the distance 6 = 5. 
The BIC energy is e^s = —0.984. (b) Vector plot of the electron 
current density in the BIC. The calculations were carried out for the 
parameter a = \/2. 


the energy Cbs and the defect potential Vbs- The electron 
density defined as |4'(a;,j/)p is shown in Fig. [^a). The 
wave function is seen to be strongly localized near the 
defect. In contrast, the edge states with very close energies 
extend along the boundary to infinity. 

An interesting property of the BIC in 2D TI is the exis¬ 
tence of an electron current. The current density can be cal¬ 
culated using the current operator for the BHZ model [I6| . 
Direct calculations show that in the BIC there is a nonzero 
current whose direction is locked to the spin similarly to the 
current in the edge states. The vector plot of the current 
density is shown in Fig. ib). 

The presence of the circulating current can be interpreted 
as an indication that the BIC is largely formed by the helical 
edge modes circulating around the defect. In this regard it 
should be noted that actually there are two BICs since the 
state we have found exists for each spin direction. In these 
states the currents are directed oppositely. However, only 
one state at a given defect can be occupied by an electron. 
Therefore around each of the defects, which forms the BIC, 
there is the electron current, the direction of which is locked 
to the spin of the trapped electron. 


4. Conclusion 

The bound states embedded in the energy continuum is a 
nontrivial quantum phenomenon which attracts increasing 
interest though they remain fragile objects in spite of the 
progress in the experiments on nanostructures. Particularly, 
very special conditions are required for their realization. In 
this paper we have shown that owing to unique properties 
of the 2D TIs it becomes possible to create the BICs in the 
potential of a quite simple form. 

The BIC can arise at a nonmagnetic defect coupled to 
the edge states in the 2D TIs. The mechanism of the BIC 
formation is due to the interference of two resonances that 
exist in the local density of the edge states flowing around 
the defect. The fact that even the short-range potential 
of the defect creates two bound states in the bulk of the 
system is a specific property of 2D TIs, which is due to the 
presence of two mechanisms of the bound-state formation. 
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In topologically trivial crystals, this property is absent. 
The resonances are formed as the defect is coupled with 
the boundary. 

To our knowledge such a simple potential inducing 
the BICs was proposed only in Ref. jT9] within a one¬ 
dimensional two-particle Hubbard model with an attrac¬ 
tive impurity potential. The emergence of the BICs on 
the surface of three-dimensional TIs was noticed in numer¬ 
ical calculations of the electron states in the corner of a 
“L”-shaped potential [20] . 

The system we have considered in this paper favorably 
differs from above systems since it is quite realizable experi¬ 
mentally. Our calculations are based on the two-band BHZ 
model that well describes the real systems. Short-range 
potential is a common minimal model of a defect. Another 
merit of this model is that it has allowed us to carry out 
explicit analytical solution of the problem. 

We have found that the BIG can arise if the distance 
between the defect and the boundary exceeds a threshold 
value. The BIG emerges at that defect with the definite 
potential the value of which depends on the distance from 
the boundary. In the future, it would be interesting to con¬ 
sider also the defect potential of larger radius which creates 
a larger number of the bound states and correspondingly 
a larger number of resonances which could facilitate the 
conditions for the BIG to appear. 

We have shown the electron density in the BIG is local¬ 
ized around the defect and has a dip in the center. We 
have also found that the electron current flows around the 
defect. The current direction is locked to the spin which 
indicates that the BIG has a helical structure. 
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Appendix A. Matrices K.{e;v,b) and jF(e;b) 

The explicit equation defining the matrices JC{e; v, b) and 
T (e; b) are as follows: 


OO 

JC{s;v,b)=J g 


_4aeAi(e,/c) 

ioo 

r dp v{k,p) 

I 2TriA{e, k,p) 


'Dq{£, k)A'{e, k)B{£, k) 


Di{e,k,p) 


, (A.l) 


= ^X>o(e,fco(e);&), (A.2) 

where — I) -I- and I?o(£, k; b) denotes the 

matrix 

g-bpi 

X>o(e, fc; b) = - Dq{s, k-pi) - Do{£, k-p2). 

Pi P2 

(A.3) 


Other matrices, Vq, 'Di,A'{£, k), and B{£, k; b), in the inte¬ 
grand are defined as follows: 


/ a22(e, fc,p) ai2{£,k,p) \ 
\-a 2 i{£,k,p) -aii{£,k,p)) ’ 


(A.4) 


f a22i£,k,p) -ai2{£,k,p)\ 

\-a2i{£,k,p) aii{£,k,p) )' 

f ai2(e, fc,P2) -ai2{£,k,pi)\ 
\-a22{£,k,p2) a22(e, fc,Pi) ) ’ 


(A.5) 

(A.6) 


B{£, k] b) = 

fv{k,pi)a22{£, -v{k,pi)ai2{£, k,pi)e~'’P^\ 

\y{k,p2)a22i£,k,p2)e~^P^ -v{k,p2)ai2{£,k,p2)e-^^^) ’ 

(A.7) 

where aij{£, k,p) are the elements of the matrix [£ — h{k,p)]: 

aii=£ + l-k‘^+p'^, ai 2 = -a{k+p), .. . 

a 2 i = —a{k — p), 022 = e — 1 + 

In addition, the scalar functions Ai(e, A:), A(e, fc,p), and 
Pi,2(,£,k) have the form: 


Ai(e, k) = a 22 {£, fc,pi)ai 2 (e, k,p 2 ) 

- ai2(e, A:,pi)a22(e, A:,P2) , (A.9) 


A{£,k,p) =-[p^-pl{£,k)\ [p^-pl{£,k)\, (A.IO) 

with 


Pi,2(e, fc) = \Jk'^ + a^/2 -l± 

and Repi_ 2 (e, k) > 0. 


4)/4-te2 

(A.ll) 
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